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Abstract. Using a ladder-rainbow kernel previously established for the soft scale of light quark 
hadrons, we explore the extension to masses and electroweak decay constants of ground state 
pseudoscalar and vector quarkonia and heavy-light mesons in the c- and b-quark regions. We make 
a systematic study of the effectiveness of a constituent mass concept as a replacement for a heavy 
quark dressed propagator. The difference between vector and axial vector current correlators is 
^> examined to estimate the four quark chiral condensate. The valence quark distributions, in the 

pion and kaon, defined in deep inelastic scattering, and measured in the Drell Yan process, are 
investigated with the same ladder-rainbow truncation of the Dyson-Schwinger and Bethe-Salpeter 
equations. 
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g DYSON-SCHWINGER EQUATIONS OF QCD 

,__! A great deal of progress in the QCD modeling of hadron physics has been achieved 

> through the use of the ladder-rainbow truncation of the Dyson-Schwinger equations 
(DSEs). The DSEs are the equations of motion of a quantum field theory. They form 

^ an infinite hierarchy of coupled integral equations for the Green's functions (n-point 

cn functions) of the theory. Bound states (mesons, baryons) appear as poles in the appro- 

lO priate Green's functions, and, e.g., the Bethe-Salpeter bound state equation appears after 

O taking residues in the DSE for the appropriate color singlet vertex. For recent reviews 

,__! on the DSEs and their use in hadron physics, see Refs. [1, 2, 3, 4]. 
L| In the Euclidean metric that we use throughout, the DSE for the dressed quark 

•^ propagator is 

X 

^ S{p)-^ = Z2i/p + Z^m{^)+Zij g^D^y{p-q)-Y^S{q)V^{q,p) , (1) 

where D^^{k) is the renormalized dressed-gluon propagator, Y^y{q,p) is the renormal- 

ized dressed quark-gluon vertex. We use / to denote / d^q/{2nY with A being the 
mass scale for translationally invariant regularization. The renormalization condition is 
S{p)^^ = iy- P + tn{ii) at a sufficiently large spacelike /i^, with m{ii) the renormalized 
mass at renormalization scale jU. We use /i = 19GeV. The Z,(/i, A) are renormalization 
constants. 



TABLE 1. DSE results [10] for pseudoscalar and vector meson masses and electroweak decay 
constants, together with experimental data [11]. Units are GeV except where indicated. Quantities 
marked by t are fitted with the indicated current quark masses and the infrared strength parameter 
of the ladder-rainbow kernel. 
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calc 


3 - 6 MeV 
5.5 MeV 


80 -130 MeV 
125 MeV 


(0.24 GeV)3 
(0.241 GeV)31' 


»% Iti mK /k \ mp fp \ m*^ f^ m^ fy 


expt 
calc 


0.138 0.131 

o.i38f o.or 


0.496 0.160 
0.497^" 0.155 


0.770 0.216 
0.742 0.207 


0.892 0.225 
0.936 0.241 


1.020 0.236 
1.072 0.259 



Bound state pole residues of the inhomogeneous Bethe-Salpeter equation (BSE) for 
the relevant vertex, yield the homogeneous BSE bound state equation 
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where K is the renormalized qq scattering kernel that is irreducible with respect to a pair 
of qq lines. Quark momenta are q^ = q + r\P and q^ = q— {\ — r\)P where the choice 
of r\ is equivalent to a definition of relative momentum q\ observables should not depend 
on r\. The meson momentum satisfies P^ = —M^. 

A viable truncation of the infinite set of DSEs should respect relevant (global) 
symmetries of QCD such as chiral symmetry, Lorentz invariance, and renormal- 
ization group invariance. For electromagnetic interactions and Goldstone bosons 
we also need to respect color singlet vector and axial vector current conserva- 
tion. The rainbow-ladder (LR) truncation achieves these ends by the replacement 

K{p,q;P) -^ -47t a^ff {k^)D^i^l^{k)YY^ ® y7v along with the replacement of the DSE 

kernel by Zig'^D^v{k)ry {q, p) -)■ An tteff (^^) ^Jlv'^(^) 7v y where k^p-q, and tteff (^^) 
is an effective running coupling. This truncation is the first term in a systematic ex- 
pansion [5, 6] of K; asymptotically, it reduces to leading-order perturbation theory. 
These two truncations are mutually consistent: together they produce color singlet 
vector and axial-vector vertices satisfying their respective Ward identities. This ensures 
that the chiral limit ground state pseudoscalar bound states are the massless Goldstone 
bosons from chiral symmetry breaking [7, 8]; and ensures electromagnetic current 
conservation [9]. 

We employ the ladder-rainbow kernel found to be successful in earlier work for light 
quarks [8, 10]. It can be written a^ffik^) = a^^{k^) + a^^{k^). The IR term implements 
the strong infrared enhancement in the region < A:^ < 1 GeV^ required for sufficient 
dynamical chiral symmetry breaking. The UV term preserves the one-loop renormaliza- 
tion group behavior of QCD: aeff(^^) —^ as{k^y^°°'P{k^) in the ultraviolet with Nf = 4 
and Aqcd =0.234 GeV. The strength of a^^ along with two quark masses are fitted 
to (qq), mj^i]^. Selected light quark meson results are displayed in Table 1. The in- 
frared kernel component is phenomenological because QCD is unsolved in such a non- 
perturbative domain. To help replace such phenomenology by specific mechanisms, it is 
necessary to first characterize its performance in new domains. 



TABLE 2. Calculated masses and electroweak decay constants for ground state pseudoscalar and 
vector heavy-light mesons, together with experimental data [1 1], all in GeV. In the rows labelled calc 
M, the heavy quark is described by a constituent mass fit to the lightest pseudoscalar (marked by t). 
In the rows labelled k^in, the heavy quark is dressed through the DSE with an infrared suppression 
of the gluon momentum as described in the text. No such infrared suppression was applied to the 
dressing of the light quark and the binding kernel. The bracketed values for Bs and B^ indicate that 
the entire a^^{k^) term of the kernel was eliminated, as in earlier work, and this corresponds to a 
^min value that is about 20% larger than the value used for B. 
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expt M 
calc M 

with kynin 


1.86 
1.85f 
1.88 


2.01 
2.04 


1.97 
1.97 
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2.11 
2.17 


5.28 

5.27+ 
5.15 
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5.41 
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expt f 
calc f 

with krain 


0.222 
0.154 
0.260 


? 
0.160 


0.294 
0.197 
0.275 


7 
0.180 


0.176 
0.105 
0.265 


7 
0.182 
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0.144 

(0.164) 


7 
0.20 
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0.210 

(0.453) 


7 
0.18 



HEAVY QUARK MESONS 



In Table 2 we display the results for the heavy-light ground state pseudoscalars and 
vectors involving a c-quark or b-quark. We use DSE solutions for the dressed light 
quarks. If a constituent mass propagator is used for the heavy quark, with the constituent 
mass obtained from a fit to the lightest pseudoscalar, the various meson masses are 
easily reproduced. The constituent masses found this way are M^°"** = 2.0 GeV for 
the c-quark, and M^°"^ = 5.3 GeV for the b-quark. To compare with what is known, 
we take the quark current masses [11] m^ = 1.2±0.2 GeV, and m^ = 4.2±0.2 GeV at 
scale /i = 2 GeV and use the quark DSE to run the masses into the timelike region 
where the meson mass shells are located. If all the meson momentum runs through 
the heavy quark, the mass function from the DSE would suggest an effective quark 
mass M^^^{p^ ~ —M^) where M is the meson mass. These DSE masses reproduce the 
previously obtained values for M^?^^ within 10%. In this sense, heavy quark dressing 
is well summarized by a constituent mass. However, the electroweak decay constants 
obtained from the constituent mass approximation are 30-50% below the experimental 
values. Moreover, within this LR model, quark dressing is not a minor effect because the 
use of fully dressed quark propagators, both heavy and light, does not yield a physical 
bound state solution for these heavy-light states involving a c-quark or b-quark. 

The results for equal quark mesons (quarkonia), displayed in Table 3, show a different 
perspective. With the same fitted constituent masses, the quarkonia masses are well 
reproduced, and again the electroweak decay constants /m are too low by some 30-70 
%. Use of dynamically dressed propagators removes almost all of this deficiency and the 
decay constants are well reproduced. This improvement provided by dynamical dressing 
of c- and b-quarks is persistent and systematic. When the dressing is progressively 
introduced into all three stages (bound state solution, normalization loop integral, and 
then the loop integral for /m), the result always increases towards the experimental value. 

Our results for decay constants /m of both heavy quarkonia and heavy-light mesons 
indicate that a constituent mass approximation, even for b-quarks, is surprizingly inad- 



TABLE 3. Calculated masses and electroweak decay constants for ground state pseu- 
doscalar and vector quarkonia, together with experimental data [11], all in GeV. Results 
employ either a constituent mass treatment of the heavy quark, labelled Mp°"^ or a dy- 
namical dressed propagator, labelled ZR^^(/?^). The last row shows the results from an 
effective smooth infrared cut-off of the gluon momentum applied to both binding from 
the BSE kernel and to quark propagator dressing within the DSE. 





Mn, 
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Mj/xf, 


fj/xif 


Mn, 


hb 


Mt 


/t 


expt 

calc with M™"' 

calcwithEg^E(-^2) 


2.98 
3.02 
3.04 


0.340 
0.239 
0.387 


3.09 
3.19 

3.24 


0.411 
0.198 
0.415 


9.4? 
9.6 
9.59 


7 
0.244 
0.692 


9.46 
9.65 
9.66 


0.708 
0.210 
0.682 


with kmim 


3.04 


0.371 


3.24 


0.398 


9.59 


0.630 


9.66 


0.621 



equate. This is strongly associated with the very weak binding of such mesons which 
magnifies the effect of a weak momentum dependence in the heavy quark mass func- 
tion M^^^{p^) and the quark field renormalization function Z^^^(p^) in the time-like 
domain where the mass shell is to be formed. A mass shell will not form unless all 
the momentum dependence conspires to push the eigenvalue of the BSE kernel to cross 
unity. 

With increasing quark mass, the quarkonia states become smaller in size and the ul- 
traviolet sector of the ladder-rainbow kernel becomes more influential. The size of the 
dressed quark quasi-particle also decreases with quark mass so that the entire heavy 
quarkonia dynamics becomes increasingly dominated by the ultraviolet sector of the 
kernel as dictated by pQCD. Although a ladder-rainbow model of ^bse has no depen- 
dence on quark or meson masses apart from quark propagators, our results for heavy 
quarkonia are quite good; one can confirm that the dynamically generated size of the 
meson wavefunction naturally concentrates the physics increasingly within the ultravio- 
let sector of the kernel. 

For heavy-light mesons, with increasing heavy quark mass, the infrared sector of the 
dynamics experienced by the light quark quasiparticle can be expected to remain dictated 
by the same soft QCD physics that determined our model kernel. The question then 
arises as to whether the self interaction of heavy quark to produce a quasiparticle, and 
its coupling via gluon exchange with the light partner, are adequately described by such 
a kernel. In a schematic analysis of ^bse beyond ladder-rainbow truncation [6, 12, 13] 
it has been shown that an important contributor to its strength of attraction is dressing 
of the quark-gluon vertex. This mechanism will decrease with increasing quark mass. 
However our ladder-rainbow effective kernel was established for light quarks and does 
not contain such a decreasing strength. 

For this reason we investigated a schematic suppression of the infrared sector 
of our kernel with increasing quark mass niq < nij,, using the form a^^{k^) = 
f{mq{^))a}^{k^) where f{mq{}i))= 1 -0.624m^(/i)/m^(/i), the renormalization 
scale of the model is /i = 19 GeV, and nij^di) =3.8 GeV. This smooth infrared sup- 
pression can be recast as an effective sharp lower cutoff k^nm for the gluon momentum 
and we estimate such a cutoff by equating integrals of the kernel: /^ a^^{k^) = 
Ik>k ■ 05^^ (fc^). The value of kj^m will be discussed below. This a^^ was used for both 
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FIGURE 1. Left Panel: The effective infrared gluon momentum cutoff versus pseudoscalar quarkonia 
mass. Right Panel: The pseudoscalar quarkonia electroweak decay constant versus quark current mass. 
Solid curve is the is from the original LR kernel without IR suppression; dashed curve is with m^- 
dependent IR suppression as described in the text. 



binding and quark dressing for the quarkonia states, and the results are shown in the 
indicated row of Table 3. Only the decay constants show a response and they decrease 
by less than 10%, as shown in Fig. 1. 

This schematic infrared suppression is not well defined for a meson with two distinct 
quark mass scales. So for initial investigation we employed the suppressed kernel for 
dressing the heavy quark and the original kernel for the dressing of the hght quark and 
the binding to form the meson. The results displayed in Table 2 show that physical D, D^ 
and B states are produced this way, whereas that was impossible with a b-quark dressed 
via the full kernel. We have not investigated the vector states this way. The pseudoscalar 
masses are quite reasonable; the D andZ)^ decay constants are closer to experiment than 
the values obtained with the constituent mass approximation. We have not tried to fine 
tune this suppression phenomenology; the /g value could be improved. The bracketed 
results for B^ and Be were produced by removal of the entire infrared component of our 
kernel, not by necessity, but for temporary convenience. 

The results indicate that our ladder-rainbow kernel, determined by essentially chiral 
quark physics, has implicit infrared strength that should be systematically reduced as 
quark mass is raised into the b-quark region. One may think of this as implicit quark- 
gluon vertex dressing in QCD that ebbs away, as it must, with increasing quark mass. 
Another viewpoint is provided by the analysis of Brodsky and Shrock [14] that, because 
QCD color confinement limits quarks and gluons to be found only within hadrons, their 
dynamics is limited by a universal maximum wavelength characteristic of hadron size. 
Thus there is a minimum momentum for such virtual fields. Since the size of hadrons 
is a dynamically generated scale, exact QCD dynamics could be expected to automat- 
ically generate a suppression of modes below a minimum momentum. However with 
truncations and models of QCD, such as the one we deal with, an infrared momentum 
cutoff may need to be a model element to represent explicitly excluded higher order 
processes. The minimum momentum must increase with mass of the hadron or dressed 
quark/gluon quasiparticle. A Compton wavelength translates to a rest mass, and thus we 



expect a minimum momentum to be proportional to mass. In Fig. 1, we display the effec- 
tive sharp kjnin as a function of quarkonia mass that corresponds to our smooth infrared 
suppression investigation. The shaded band estimates the associated uncertainties. The 
increase in proportion to the mass is as expected. 

LEADING NPQCD SCALE AND THE FOUR QUARK 

CONDENSATE 

Quark helicity and chirality in QCD are increasingly good quantum numbers at short 
distances or at momentum scales significantly larger than any mass scale. One manifes- 
tation of this is that, for chiral quarks, the correlator of a pair of vector currents is iden- 
tical to the corresponding correlator of a pair of axial vector currents to all finite orders 
of pQCD. Non-perturbatively, the difference of such correlators measures chirality flips, 
and the leading non-zero ultraviolet contribution identifies the leading non-perturbative 
phenomenon in QCD. This is the four quark condensate [15]. The LR kernel can produce 
the difference correlator as a vacuum polarization integral in momentum space, where 
the propagators are dressed and the vector and axial vector vertices are generated in a 
way consistent with the symmetries. We use the large spacelike momentum dependence 
to extract the leading coefficient or condensate. 

The vector current-current correlator is formulated as the loop integral 



n%{P) = Jd'xe'''-'{0\Tj^{x)j+iOm = - j'^Tr{y^S{q+)Yl{q,P)S{q^)} , (3) 



A 

V 



where A indicates regularization, e.g., by the Pauli-Villars method, and T\, is 
the dressed vector vertex. The axial vector correlator is formulated in an anal- 
ogous way and we directly calculate the difference correlator which does not 
require ultraviolet regularization. With nKy(P) = (P^S^v — ^/i^v)n^(i'^), and 

I\^^{P) = {P'^5^v-P^Pv)Yl'^{P^)+PpiPv^t{P^), the quantity of interest here is 

The leading non-perturbative contribution to D^ starts with dimension d = 6 and 
involves the four-quark condensate in the form [16, 17] 

Our numerical calculation of P^Y\^^^[P^) identifies a leading ultraviolet constant rea- 
sonably well. The four quark condensate {qqqq) extracted via Eq. (4) is 65% greater 
than the common vacuum saturation assumption {qq)^ at the renormalization scale 
/i = 19 GeV used in this work. The low P^ limit provides a reasonable account of the 
first Weinberg sum rule [18, 19]: P^Yl'^j~^{P^)\p2^Q = -f^ in the /;, = 0.0924 GeV 
convention. This limit is due to Ilj only and we obtain f^ — 0.09 GeV that way. Our 
results are consistent with the second Weinberg sum rule [18] P^Tlj^^{P^)\p2^^ — 0. 
The Das-Guralnik-Mathur-Low- Young sum rule [20] relates Jq dP'^ P'^U\-^'^{P'^) to the 



strong component of mj^± — m^o. We obtain 4.86 MeV for this mass difference in com- 
parison with 4.43 ±0.03 from experiment. The chirahty-flip ratio n^"'^(p2)/n^+^(/'2) 
identifies a scale of ~ 0.5 fm for the onset of non-perturbative dynamics [21]. This DSE 
model in LR format produces the same scale as obtained in both a lattice-QCD calcula- 
tion and the Instanton Liquid Model [22]. 



VALENCE QUARK DISTRIBUTIONS IN THE PION AND KAON 

Data for the momentum-fraction probability distributions of quarks and gluons in the 
pion have primarily been inferred from Drell-Yan [23, 24, 25] and direct photon produc- 
tion [26] in pion-nucleon and pion-nucleus collisions, and semi-inclusive ep — )■ eNX re- 
actions [27]. For a recent review of nucleon and pion parton distributions see Ref. [28]. 
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FIGURE 2. Left Panel: Pion valence quai'k distribution evolved to (5.2 GeV)^. Solid line is the full DSE- 
BSA calculation [29]; dot-dashed line is the semi-phenomenological DSE-based calculation of Hecht et 
al. [30]; experimental data points are from [25] at scale (4.05 GeV)^; the dashed line is the recent NLO 
re-analysis of the experimental data [31]. Right Panel: The ratio of u-quark distributions in the kaon and 
pion. The solid line is our preliminary result from DSE-BSE calculations [32, 29, 28]; the experimental 
data is from [33, 23]. 

Lattice-QCD is restricted to low moments of the distributions, not the distributions 
themselves [34]. IVIodel calculations of deep inelastic scattering (DIS) parton distribution 
functions are challenging because it is necessary to have perturbative QCD features 
(including the evolution of scale) coexisting with a covariant nonperturbative model 
made necessary by the bound state nature of the target. Aspects of chiral symmetry have 
led to DIS calculations within the Nambu-Jona-Lasinio model [35, 36, 37] but there are 
a number of difficulties with this approach [36, 37], among them a point structure for 
the pion BS amplitude at low model scale and a marked sensitivity to the regularization 
procedure due to the lack of renormalizability. Constituent quark models have also 
been employed [38, 39], with the difficulties encountered in such studies considered 
in Ref. [40]. Instanton-liquid models [41] have also been used. In these approaches, it 
is difficult to have pQCD elements join smoothly with nonperturbative aspects. All of 



these issues can in principle be addressed if parton distribution functions can be obtained 
from a model based on the DSEs, and in particular styled after the approach we have 
previously made for the related pion electromagnetic form factor [42, 43]. 

In the Bjorken kinematic limit, DIS selects the most singular behavior of a correla- 
tor of quark fields of the target with light-like distance separation z^ ~ 0. With inci- 
dent photon momentum along the negative 3 -axis, the kinematics selects z"*" ~ z^ ~ 
leaving z^ as the finite distance conjugate to quark momentum component xP^ , where 
X = Q^/2P -q is the Bjorken variable, q^ = —Q^ is the spacelike virtuality of the pho- 
ton, and P is the target momentum. To leading order in the operator product expansion, 
the associated probability amplitude qf{x), characteristic of the target, is given by the 
correlator [44, 45] 

qfix) = ^Jdz-e^^'^^-{n{P)\xfrf{z-)r+Wfm^{P)) , (5) 

where / is a flavor label. Here the 4-vector components that arise naturally are 
a± = ((3*^±a^)/\/2. This probability amplitude is invariant, it can easily be given a 
manifestly covariant formulation, and its interpretation is perhaps simplest in the infi- 
nite momentum frame where qf{x) is the quantum mechanical probability that a single 
parton has momentum fraction x [46]. Note that qf{x) = —qA—x), and that the valence 

quark amplitude is q^Ax) = qf{x) —qAx). It follows from Eq. (5) that /q dxq^Ux) = 
{7r{P) \jt {0)\7r{P)) /2P^ = Fj^(0) = 1. Approximate treatments should at least preserve 
vector current conservation to automatically obtain the correct normalization for valence 
quark number. 

In a momentum representation, qf{x) can be written as the evaluation of a special 
Feynman diagram [44, 47] 

^M = \[ 5ik-^-xP+)\.v,A[y+G{k,P)] , (6) 

where trcd denotes a color and Dirac trace, and G{k,P) represents the forward ^-target 
scattering amplitude. In ladder-rainbow truncation, which treats only the valence qq 
structure of the pion, Eq. (6) yields the explicit form 



G}(x) = y tv,i[r^{p,p)s{p)r+{p;x)sip)r^{p,p)sip-p)] , (V) 



where r+(p;jc) is a generalization of the dressed vertex that a zero momentum photon 
has with a quark. It satisfies the usual inhomogeneous BSE integral equation (here with 
a LR kernel) except that the inhomogeneous term is 7+ 5(/?+ — xP+). This selection of 
LR dynamics exactly parallels the symmetry-preserving dynamics of the corresponding 
treatment of the pion charge form factor at ^^ = wherein the vector current is con- 
served by use of ladder dynamics at all three vertices and rainbow dynamics for all 3 
quark propagators [42, 43]. Here the number of valence M-quarks (or d) in the pion is 
automatically unity since the structure of Eq. (7), along with the canonical normalization 
of the qq BS amplitude Y^{p,P), ensures /q dxq\^{x) = 1 because /q dxT^{p\x) gives 
the Ward Identity vertex. 



Eq. (7) is in Minkowski metric so as to satisfy the constraint on p+, but LR dynamical 
information on the various non-perturbative elements such as S{p) and Yti{p,P) is 
available only in Euclidean metric [10]. Since qf{x) is obtained from the hadron tensor 
W^^ which in turn can be formulated from the discontinuity T^^{e) —T^'^{—e), we 
observe that all enclosed singularities from the difference of Wick rotations cancel 
except for the cut that defines the object of interest. With use of numerical solutions 
for dressed propagators and BS amplitudes, that give an accurate account of light 
quark hadrons, our DIS calculations significantly extend the exploratory study made in 
Ref. [30]. That work employed phenomenological parameterizations of these elements. 

In Fig. 2 we display our DSE result for the valence w-quark distribution evolved to 
2^ = (5.2 GeV)^ in comparison with kN Drell-Yan data [25] with a scale quoted as 
Q^ > (4.05 GeV)^. We also compare with a recent NLO reanalysis of the data at scale 
Q^ = (5.2 GeV)^. The distribution at the model scale Qq is evolved higher by leading 
order DGLAR The model scale is found to be Qq = 0.57 GeV by matching the x" 
moments for n = 1 , 2, 3 to the experimental values given independently at (2 GeV)^ [48] . 
Our momentum sum rule result /q dxx{uji + dji) = 0.74 at Qq clearly show that in a 
covariant approach the retardation effects of one gluon exchange assign some of the 
momentum to gluons. The corresponding momentum sum for the kaon is 0.76. 

The ratio uk/u^ measures the dynamical effect of the local environment. In the kaon, 
the M-quark is partnered with a significantly heavier partner than in the pion and this 
shifts the probability to relatively lower x in the kaon. Our preliminary DSE model cal- 
culation [32, 29, 28] is shown in Fig. 2 along with available Drell Yan data [33, 23]. 
Here we include only the leading two invariants of the pion BS amplitude, E{q^P) 
and F{q^P), where q is qq relative momentum. For both amplitudes only the lowest 
Chebychev moment in ^ ■ P is employed . This variable does not occur in static quan- 
tum mechanics, nor in the Nambu-Jona-Lasinio point-coupling field theory model [35] 
which also neglects the q^ dependence. We do not make such a point meson approxima- 
tion here; the q^ dependence comes from the BSE solutions. Nevertheless, the essential 
features of the ratio uk/uti are adequately reproduced by a generalized Nambu-Jona- 
Lasinio model [28]. 
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